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Awe HISTORY OF GEOMAGHETICS 


The beginnings of the study of geomagnetism lie back in 


prehistory when magnetic attraction between iron and certain 


Minerals was first observed. Exactly when “his phenomenon 
was first noticed is not known, but the proverties of magne- 
mee, -hen calied lodestone, appeared in Greek literature 


around 600 B. C. (Bre¢nnan and Davis). [Ref. 1] 
Snaonan [| Ret. 2] 2ndicates that the directional property 

of magnets was known and ased in Burop2 prior to 

and possibly in China before then. EF. N. Pa 

Mee-a "Te 35 an interesting fact that the ancient walis of 


Peking were lined up with a 


a a 
Geographic north, a difference at that time of about 19. We 
may presume that the surveyor found it easier +t0 work with 
his compass needle by day than tc sight on the pole ster by 
Cee SGU baste Ss Jeotsd- 


Megat. «6S property allows the uss 
se 


Pe -1¢ f2eld tor navigational purpo 


Beene M@od—f2tceentn can =ury it was determined ir Europe 
ya. the Magnetic compass does not point to true north. Thea 


eg] be Weer true no=zth and the direction indicated bv «he 
compass is now known as magnetic declination by the geophy- 
Stcist and as variation by the navigatecr. 

nema gie:.ic £~eld dip, Sr magqne*=ic inclination is th 
angle, in a vertical plane, between the horizontal and “<the 
Mecect.on Of ‘the earth's Magnetic field vector. i+ We 
Observed in 1544 by an instrument maker in Nuramberg named 
Hetumenn, and again hy Robert Worman if London in 15861 

These discoveries or observations gave rise + 


te) 
of geéomacnetics as a specialty which gained its cornersto: 


+? 
t- 





with William Gilbert in 16009. Ai~e=5 COMmpareng Lic =xperi- 
mental results with the previous work of others such as 
merMan, Gilbert, in his book, "De Hagnete,™" concluded «hat 
"Magnus Magnes jpse est globus terrestris (*he earth globe 


itself is a great magnat)" [Ref. 2]. It Sse0his copreceépt, 


emt the earth is »tself a magnet, ‘hat is the basis of ths 


Science of ape. CS. 


Gilbert felt tha the Barth's maqnetism must reéenain 
constant except for Ns cess changes, but i+ was soon 
determined that this was not the cas¢. A 'secular varia- 
meom Of the Earth'< was found co exist. 


Shorter term changes in the azesmagnetic field were 


observed and it was eventually realized that geomagnetism is 
dynamic. In 1722 George Srtaham discovered «that daily, of 
G@eucnal, Variations exist [Ref. 2]. Bieing ~ve. ea =i 


Nineteenth century magnetic observatories began to be astab- 
iene d to record «the changes in the geomagnetist field in 2 


systematic fashion (Knecht) [Ref. 4]. 


Be EARTH'S MAGNETIC FIELD 


There are various ways of breaking down the consti- 
hg parts of the geomagnetic field. One way is <9 divide 
the field 


Deing this yields these three parts: internal, crustal, and 


erties ocr: darstenee Fron the center st the earth. 


rw 
se) 


exeeroe!l (RPGL) (Bef. 5). Tiewamtsenalerteld ofFicinaetes in 
m@emce5e =region and if the more stable field, conzainiag 
only extremely low frequency temporal var 
crustal (or anomalous) fi2ld arises from medificazions nad 
a Uneomeeneeaneal fdscld by aeterials aed structures in the 
fac th "Ss “@F wst . These Vatlations ape net spa tially consfan 


and give rise to some of what is known as gqeclogical 





variations. The external field is “ke ost dynamic and 


arises from many sources including + ction between 


% 


n a 
the solar wind and the terrestrial magnetic field. 
ie 


Another way of d2scribing 


is @ | 
O 
i} 
iD 
ty 
ct 
Ww 
O 
th 
ct 
rye 
(D 


on 
geomagnetic field is by time variatior. This division is 


accomplished by considering that opart of the field which 


steady rie¢id and what is left as the variation Ffisld 


(Knecht). f[Ref. 4] 


The st2zady field consists of the above named 
internel field, alse referred to as the main Field. Slow 


td 
(tT 
(D 
=| 
rey) 
! 

| 


Variations i field with periods cf years or longer 


are referred to as secular variations. 


Sa ea ae oe ae ——P Ep ae GE a 


2. Wodeils of the Main Fiela 





Various nodels of the main geomagnetic Field make 
use of a geocentric dipole. Gauss, in 1839, den a 
that, asa fairly good first approximation to the geomag- 
netic field, the field of a uniformly magnetized sp ( 
points outside the sphere) is equiv2zlent to the field o 
magnetic dipole located at the Gemeor wort” =re enhers 
(Jacobs). {Ref. 6] 

The simplest cf the present approxinations of the 
@eeMagnetic Fieid ws that of a short bart magne: or dipol 
located at the center cf the earth with an axis inclin 
eeprOoxm@m@ately 711.5° from the Barth*s axis of rotation. The 
sense of the field lines is from south to north (Figure 
eek ice 


ctr 


Mie exis of this dipole in®fersects the sarth at the 
S€enagre=2c north polé, 738.59N, 291.09EFE im geographic ccor- 


865°S, 11080%. 


~ 


dinates, and at the geomagnetic south pole, 
The moment of the geomagnetic dipol= is 9.1 


c 
It is these poles that area used to define the aqecmacnetic 


3 
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fe 
SS , ee ee a 
Figure 1.1 Dipole Appearance of Geomagnetic Field 
coordinate system (Knecht) [Ref. 4]. The geomaqneztic coor- 


dinate system is a spherical polar system sinilar> *o «he 

geographic coordinat2 system with 2 geéomagrstic equator 

defined 90 degrees away from ¢€ith2r geomagnetic role in 
Cc 


latitude. This tilted geocentric dipole model describes the 
Main geomagnetic field to an accuracy cf about 10%. 

In 1940 Chapman and S3eartals dafined in of f-centsr 
dipole in the earth's interior, c 
Mies edipeto “sso displaced 0.0695 
Magnitude from the center and in the 
(oon IDOI OE a(igeegneaphic coordina 





hasweeamtets=ctions Of the eccentric dipole axis at the 
Semen s Suerace are 81.99N, 84.79W and 75.9°9S, 120.4908 
(Fiaqure 1.2) (Haymes) [Ref. 8]. This approximation is accu- 


rate to within a few percent. 






Geomagnetic 
equator - 
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Pigure 1.2 Eccentric Dipole Model of Geomagnetic Field 


The field has additionally been modelsi *o an accu- 


Eeemwet abou. 1 by detezhinining Seussian Gostficients by a 


(D 


tn 


rt 


least-squares fit of experimental measurements o 





geomagnetic field. These coefficients 
ical harmonic series representing the sc 
the field. This accuracy implies that th 
bution +o the total main field is at leas 
oO 

The International Geomagnetic Reference Field (IGRF) 


Meetas vVelues which difitsc by only parts oper thousand fron 


measured values. 


3. Sources of the Seonagnetic Field 


[= a ee Se =o Ea aE aa 


There are various elements thatc u 
geomagnetic field, some external to the earth's s 
some internal. AS previously mentionne 
contributions make uv only a small fra 

teld, playing a more important role in the variation field 
Thes2 external sources include current systems in 
the earth's upper atmosphere affected by sola Ste 

netic radiation and gravitation, solar corpuscular ra 

Se, =ne interaction of solar plasma with the main fiel 
the effect of the solar interplanstary field. ‘Ref. 47 

Various magnetic surveys of the world, _neiud sag 
those ccrnducted at ground level, by airborne instruments, 
ema by seatellits,have pointed to the fact chat the largest 
source of thé ecm onemre cle ld 2S — 2r~sronal te it. 
While there exists resiiual pernmanen m 
Paeee se chusce, -his cannot be the principa 
Grec=ne]e geomagnetic field due to ¢ 
properties knewrn t0 exist in the earth's interior (Néegata 
and Ozima). {[{Ref. 9] 

Fermanent magnetism is generated by micrescooic 
electric currents, since a changing ele 
generate a magnetic fieli. Anothe 


Taeeeatcettmeid is by the Qocticn of electric charges ih 


ja 





PecTOscopic current. 


conductin 


SeEO01LC Cu 


g fluid core 


rrent systen, 


Convective motion cf 
of *he earth, 


1s considerad 


eee. ema cally 
Besu lian. of 


to he the 


source of the main field. 

The most promising present theory of the generation 
of the geomagnetic field is that of some sort of a selfi-ex- 
citing dynamo systen,. This means that the motion of 4a 
Genductor, such as the molten iron in the earwh's core, in a 
Gegretic field produces a current which in turn induces a 
Magnetic ficsld in support of the original magnetic field 
tBer. 4, 7). A very Simple model of such a dynamo is shown 
Perea gure 1.3. 





Axis of rotation 


Nonrotsting coil 
Rotating disc 
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Figure 1.3 Simple Disk Dynamo 
The original excitation er seed fieid nay be due to 
eee Oxcernali fiela@ line, perhaps from the solar intjerplane- 
Dorvemcamettc fCield. This original, poloidal, field line is 
Wemma up dae to the dirt@rential rotation (rotation rot 








Comstean« Waeeh Latituee) in the molten corse. The wound up 
line becomes an intensa2 azimuthal field which is carried 
cutward by the upwelling associated with convection and 
jmested by he action cf the Coriolis force. The twisting 
generates a helical toroidal field which, by cutward diffu- 
Sion, senerates the externally observed quasidipole gtomag- 


netic field. 
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The cemorthactOn sand interaction of *wo (se 
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1.4) or more disk dynamos carn als>d can also e¢xvlain the 
reversal of the geomagnetic field [Ref. 7, 9]. 
as 4 
| 
| ! 
| 
| | 
| 
| 
| 
| | 
| | 
| 
| 
| 
| 
| 
| 
7 | 
| [Ref g] | 
a _ 


Figure 1.4 Twin Disk Dynamo 


Regional ancmalies tha* are nondipolar (do ror 
Gen [Opn to wae dipole field) poss@bly arse freon eddy circue 


dations in the outer core { Ref. 4]. 





The magnetosphere can be defined to be that region 
(see Figure 1.5) occupied by the geomagnetic field above the 
ionosphere, a region where the field stongly influences the 
dynamics of ionized gas and charged particles (Kern) 
{ Ref. 10 ]. 

If the space surrounding ths earth were a perfect 
vacuum, the earth's magnetosphere and magnetic fieid nigqh+ 
be more or less symmetric and extend outward until it mezrged 
with, and its strength became insignificant compared <to the 
solar ard other planetary magnetic fields. Na Scey =e 
not to be the case. 

Instead of being an island in 2 perfect vacuum, «he 
earth encounters a continuous flow of hot, highly conductive 
ionized aas, or plasma, streaming outward from the sun. 
This continous stream of charged particles is cali2d the 
Solar Wind. The density of this ‘wind! near earth is on the 
order of 10 particles per cm3, and has a velocity av2rading 
300-500 km per second (Jacobs). [{Ref. 6] 


Bee tmeae Sele=z Wind end he geomagnetic fi 
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pressure. The hot plasma of the 


ad ¢xer" 
against 


Lar wind pushes 
da 
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the geomagnetic field deforming the field. A 
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solar wind greatly exceeds that of the qeomagnetic a 
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the geomagqneti field will be swept along wit 
Boome £5 10 earth radii inward, *he geomagne*ic field 
pressure will predominate, ¢xciuding t=he solar wind, this 
being the region of the maqnetosphere. 

In the intermediate region, she magnitude of the 
solar wind and geomagnetic field pressures are comparable 
and the solar wind is compressed and flows around the 
geomagnetic field. This occurs when the magnetic en=rgy 


density ahead of the plasma equals the kinetic energy 
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density of the streaming plasma. The solar wind 
meecnas point and forced to flow around the magznetcsnoh 
This region where the magnetospher= star*~s is call 


magnetopause. [Ref. 8] 


The velocity of the undisturbed solar wind is anale- 
gous to a ‘supersonic’ velociy. Pics. a SHOCK “EFont 2s 
formed between the magnetopause and the solar wind. The 


Maqnetosheath is the region of severe turbulence that exists 
between the shock front and the magnetopause. 
Since the solar wind always travels outw 


a 
sun, the effect of the wind on the earth's main field is nox 


completely symmetric, although it is almost synmetric abour 
an axis through the earth and sun. A geomagnetic tail is 


formed where the wind sweevs the geomagnetic field aiong 
Ween 2c Ch the nightside of the earth [ Ref. 4]. E=qu22 1.5 
pictorially represents the effects of the solar wind on the 


geomagnetic field. 


5. Time Variatiors of the Geomagnetic Field 


The geomagnetic field changes with tine. AS pDrevi- 
ously mentioned, very slow variations in the nain Tield with 
periods on the order cf years t9 thousands of years are 
referred to as secular variations. Secular variations ars 
Geologic or 'paleomagnetic’ in origin. Secular variations 
are not caused by a strength or orientation chand= of the 
center dipole. Paleomagnetic studies are used *o de~ermine 
mene secular “vVartLation. Geolegzc Structures, especially 


conductivity structure, may partially mask the secular vari- 
ation atone point on the earth 2s compared to that ar 


another point. 


Other time variations of the field can be cartéegor- 
=aec yimeo Giaet Yatdatzon fields and disturbed variation 
fie 1 d'se Disturbed variation fields include geomagnetic 
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ae Quiet Variation Pields 


Quiet vériation fields are thass 


due to disturbances in the interp 
2 


which vary slowly and regularly [Ref. 4). 
There are several contributing 

an2et Variation. These include the Solar 

Variation (Sq), the Lunar Daily Variation (L), 


Variation due to magnetosnheric effects. 
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MipasOlas Quset (Sq) wariacdon is hame given 


the 
to the pattern of diurnal field variation with respect t9 
=e 


solar local time which is caused by currents Wing in «he 
ionosphere (Matsushita) [{Ref. 11]. The major portion (about 
two-thirds) of the Sq field is due t9 what is referred as an 
atmospheric dynamo. High speed tidal winds are generated by 
solar heating causing convection of the upver atmospnhers 
(Ref. 4}. These winds produce a stationary current svsten 


by moving the conducting particles 2f the upper atmosnvhers 


W 


across the geomagnetic field lines. thie Gag ly Variation. i 
caused by the earth rotating under ths current systen. Th 


D 


remaining third of the Sq variation is caused by cutrents in 
the earth induced by the primary currents in the ionosphere. 
The Sq field can be shown to be latitude depéen- 


dent reaching a maximun 2t the magnetic equator where 432 


concentration of current, the equatorial electrojet, exist 
fRef. 4]. The maxXimum horizontal component intensity is 
geome 100 nT ae the equator with 25 *9 50 nT gore iikely ac 


higher latitudes. 
Pengitud:nraly, seasonal, and seclar cycle dep-n- 


dencies also occur for the Sq field. 


An ¢xanple of the quiet-day Vig SS Ot eee 
Monterey, California is shown in Figure 1.6 and is summa- 
Bigedq in Table I. This data was taken using a2 C2sium Vapor 


total field magnetometer in February, 1979. 


Toe Limes Daily Variation, L, iS 2pproxinetely 


Gtlemmer.h the*magnitude of the Sq field and exhibits 3 
Semaeaiurnal pbeliayior in luner time (Ref. 4}. Lhe Na joe 


difference is that the winds are caused by lunar-solar 


gravitational tides. The L field is dsperndent cn seasonal 
influences, lunar phase, the solar cycle, and lati«uds. 
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YABLE I 


Quiet Day Variation in La Mesa Village, February, 26, 1979 
POSITION 36°36'N, 121°51'W CLA MESA VILLAGE) 
CALCULATED VALUES USING THE 1975 U.S. CHART MODEL (WORLD DATA 
CENTER A, BOULDER, COLORADO) FOR FEBRUARY 1979: 
fe) I HCNT) ZCNT) FCNT) 
WALUE (FEB. 1979): 15.96°  60.74° 24,650 44,007 50,441 
YEARLY CHANGE aaa 5 =a, ° -4.3 -32.) -30.1 
MEASURED VALUES OF TOTAL FIELD INTENSITY FOR FEB. 25, 1979 
DATE TIME FCNT) 
2/25/79 00:44 $0,405.24 
91:48 $0,405.79 
02:52 $0,404.80 
03:56 $0,405.48 
05:00 50,408.19 
06:04 $0,409.15 
07:08 50,409.66 
08:12 ; $0,399.69 
09:16 $0,381.59 
10:20 50,363.88 
11:24 50,345.03 
12:28 590,338.52 MIO-OAY LOW 
Ls. 3< 50,356.48 
14336 50,378.86 
1§:340 50,386.98 
16:44 50,384.30 
13:51 50,400.81 
19:55 $0,396.91 
20:59 $0,401.35 
22:03 $0,400.30 
23:07 50,400.54 
2/26/79 00:11 50,406.51 
AVERAGE $0,390.50 + 21.43 (NT) 
fie ak cle Cages: = 5O,5-NT OR .24 
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A diurnal effect due to the dayside-nightside 
difference in compression by the solar wind of the geomag- 
netic field causes a small variation of the order of 3 aT 
fRef. 4}. 


bo Disturbed Variation Fielis 


veriations that appear to be the rasult of interplaneta 
environmental changes and do not possess a simp 
mei ty.. These variations include ionospheric dis 


Cc 
the aurora, geomagnetic Sstoras, and geomagnetic 


bemces produced by X-rays emitted from the solar flare. 
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yoically a few minutes in 
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SFE*s usually have a rapi 
G@iedeics, Loliowed by a slower return to normal. The antize 
See lasts On the order cf an hour (Reid). [Ref. 12] 
Auroras ar2= caused by *he precipitation of 
charged particles down magnetic fieli lines into the atmes- 


phere and can be one of the brightest visual ohenomena in 


the sky. Tivew QWemve matense and ac®ive auroras oceur wit! 
geomagnetic disturbances and grsatly increass icnizéticn as 


well as creating che spectactular visual displays. (Ref. 4] 
Geomaqnetic storms are due *o @ change in «he 

dynamic pressure of the solar wind. RP cyplcal stoeh Meqains 

Wee: a COmpressicn of the magnetic field by an increase in 

solar wind dynamic pressur2 called as 

(SCH » which increases the ie 

"gradual storm" begins with a graiual increase in fieli 

strength). The increase in field s 


of several tens of nanoTesilas (nT) and «a 
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Sex Wine tes to cases, if 2 disturbance starts with an SC but 
lacks the succeeding stages of a storm it is ref 
a sudden impulse (SI). Pollowing the SC, the fi 
compressed fer two to eight hours in the initia 
the storm. The main phase follows the initial phase. Over 
a period of hours to a day a westward ring currant is set up 
at a distance of several earth radii whos magnetic fieli 
leads to a decrease in field strength on the order of 100 
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ial. This decrease ove 
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ershoots the egiailibrium fleld strength 
and leads into the recovery phase of a dav or longe 


the field returns towards its prestorm strength as the 
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d 
currents gradually dissipate Ret. ie (Mazsushita) 
fRefs 13]. A magnitude-time graph o£ a typical gecmagnetic 


Seem 1S shown in Figure 1.7. 
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Figure 1.7 Typical Midlatitud2 Geomagnetic Stora 
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Geomagnetic micropalsations are rapid fluctua- 


tions in the surface magnetic fisli with perisds of about 


0.2 seconds to 10 minutes (frequencies about 90.0016 to 5.0 
fz} These are observed as a type of geomagnetic distur- 
bance by ground based magnetometers (Ref. 4}. 


Micropulsations will be discussed in depth later. 


6. Elements of the 
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The geomagnetic field vector is measured or charac- 
terized at any point by its direction and magnitude. This 
can be done in terms of some set of three independent paran- 
eters such as two direction angles and the magnituds, or 
three perpendicular components. [Ref. 4] 

The system of coordinates commonly 2nr 
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describing the geomagnetic field on the surface of the earth 
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is shown in Figure 1.8. 
locai (gecdetic) coordinat3s with respect *o True North. 


The varicus coordinates are referred *9 as maqnetic 
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aments and are defined as follows: 





Be Tomar Field Pntensity (the symbol F is 
also used) 


M@  Horizon=21 Component 


Mew N@gwehwascd, OF Nerth-South Componen= 

Ys EaStward, or East-West Component 

Z; Downward, or Vertical Compon:3n* 

D: Decanats2on OL MNagnem2c  fanptation. 
This is the angle between X and d and is 
measured positive eastward. 


ieee tnemematLOon Of (Gip angle. This is <he 
angie between H and B andis neasured 


positive downward. 
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II. INTRODUCTION TO MAGNETIC ANOMALY DETECTION (MAD) 


A. DEFINITION OF A MAGNETIC ANOMALY 


A Magn¢etic anomaiy is defined as any spatial variation 
Or disturbance in the geomagnetic field which is due to 
local causé¢s. Anomalies can be caused by waves, ore 
deposits, s2a nounts, and magnetized objects such as surfaces 
Ships and submarines (Anderson) [Ref. 14}. For the purposes 


of this research, magnetic anomalies due to a simple dipols 
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eld, such as those generated by subdmarines, wil De 


th 
de 


regazded as signals, while other anomalies will be regarded 


as noise and will be discussed later. 
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eee mao sORTCAL DEVELOTEENE OF THE MAGNETIC DETECTION O 
SUBMARINES 


ie wiemely Deteceion Systems 


Attempts at finding submerged submarines by sensin 
disturbances in the geomagnetic fieid date back as least as 
far as World War I. im 1918, “sarasst IMerarct, 2: the Naval 
see | 


experimental Station, New London studied the use a 
ving boats and airplanes 


Bees type Of detector f95 ase in mo 
(Sizchter). [Ret 15 ] 

MAD, Originally known as Magne+ic Airborne Detector, 

a@nteoow Called Magnetic Anomaity Detection, hai its beginr- 

1940 and early 1941. 


eld caused ov <the 


hnings as an airborne ASW sansor in late 
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Anomalies in the geomagnetic £ 


presence of submarines, are on the order of one to a few nT 


(Gaammas) in magnitude. This is quit? small compared tec “he 
Magnitude of the field itselst (30000 - 60000 xT). Two 
methods ars generally e2mployed to measure this smail 
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disturbance. One is to us? a gradisometer which measures the 
Spatial rate of change of the magnetic field or its 
@uadient. The second method, and the one presently used bv 
U.S. Navy aircraft is to use a magnetometer to directly 
measure changes in the magnitude of the magnetic field. 

The British investigated the use cf a gradiomneterz 
fer submarine detection, and by early 1941, had developed a 
two-ccil gradiometer system which could detect a submarine 
at a range of 200 feer under favorable conditions. aes 
renge was considered to b2 too small +o be of operational 
value and work on such 2 system was terminated when &@ nagne- 


tometer system showed promise (Coleman). (Ref. 16] 


By iate 1940, Lp icegeys ae Vacquier of «ne Gulf 
reseatch and Develcpment Company had devselooed a Sensitive 
saturable core magnetometer mntend=d for geophysical 
(mineral) prospecting. Th? Vacquier magnetometer became the 
basis for further MAD development. The Airborne Instrument 


Laboratcry of Columbia University continued the investiga- 
tion of means of localizing submerged submarines by MAD. 

The simplest saturable core or fluxgat> nagnetomneter 
consists of 2 saturable or ferromagnatic core around which 3a 
coil of wire is wrapped (s2¢ Figure 2.1). This coil carries 


Sescemusotdai current, f(t) Which is large endsugh *9 Ssacu- 
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eack cycle. Toe tnduccaniem o 
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Beeenche core during part o 
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the coil wili chance as a function of the magnetization o 
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the core. mye cOGe Badn=t ization, =m Lure, depends or th 
PieomatmeensOuS Current in the coil, and, if present, “he 
@xt<stral magnetic field. 

In the absence of an externally apolied magnetic 
ford senagnec Zeateon me 32 <uncesion of @&cliting current is 
symmetric around I=0 (Pigure 2.2). an e¢x*+ernal o 
2eld partaliel to the cora's axis will change the magn 


Pu mOn the cone and shift the WMagqnetization curve. [h2s 


30 














_ > _ 
| 
| 

| | 

| Excitation \ 
| Secondary | 

| 
| | 
| 
| | 
| 
| Ambient Field | 
| | 
| | 
| | 
| 
| 
| | 
| 
| 

{| 

— Se eee Se eee meneame eet ene ate ane me ne | 
Pigure 2.1 Element of Fluxgate Magnetometer 

Chang will cause an asymmetc which can b= sensed by 

analyzing the harmonic content of tha signal. The co"! cas 
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cetector alement (coil) alicned with the qeomaegne 

This was the basis of the MAD magnetometers used in 
Moe wee: EE Such as “ha AN/ASQO-1, ASO=n, aba ASO=—2 
[ REL. 16]. 


w Oo 


The fFluxgate magnetometer measures only the cempdo- 

Ment Of the excernai field parallel to the axis of the 

ferromagnetic core. hReeon cer to measure the total field “rn 
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wets fashion Lt is necess <5 cise, the farrcmagqnsic cazs 
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Figure 2.2 Magnetization Versus Exciting Current 


perpendicular (orthogonal) cores. Precision requirements 
made this type of system difficult to sresalize during World 
War II but it has since found many years of operational 
usage and has found widespreéed use ln geophysical explora- 


tion work and satellite mapping of the geomagnetic field. 
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Duzing World War II MAD provided a passive methed of 
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detection and tracking of submerged submarines. 
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surprise to the attacking aircrew. Until the 


ac 
was delivered the crew of the submerg2d submarine might no* 
even be ewarte of the aircraft hunting it. 

Operational deployment of MAD began ona limited 
basis in December, 1941 with the installation of she early 
oa s few ino@ blimp at Naval Air Station Lakehurst, New 


Jersey [Ref. 16}. 





By the end of 1942, MAD was operational on board PBY 
Cataline aircraft, nicknamed "Madcats." Even though MAD was 
Operational at this time, it was not until February, 1944 
[meeeean amctial contact Dy MAD led to the sinking of a 
submarine. 

During February, 1944, the Madcats of Patrol 

Squadron 63 were assigned to fly a MAD barrier patrol of «the 
Seees=ts Of Gibraltar. On 24 February 1944, a Madcat of 
VP-63 detected U-761 by use of MAD and commenced tracking to 
confirm the contact as 2 moving target An attack was 
conducted in Conjunction with nother Catalina, = wo 
destroyers, and eventually two other aircraft, and U-761 was 
sunk. U-392 was sunk after a similar MAD contact on 16 
March, and on 5 May, the third successful attack on a U-boat 
resulting from an initial MAD contact took place when U-731 
waS sunk (OEG No. 51, Price). (Ref. 17, 18] 
. The attacks at the Straits of Gibraltar denonstrated 
that though MAD has a limited search rate, thera ar: 
scenarios where it can be employed effectively as a search 
Sensor, such as providing a blockade across a4 restricted 
area without the presence = SUmeace, G=aare. (2EG No. 54) 
iene 19). 
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The magnetometer system in curren* operational use 
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is the optically pumped magnetometer. he optically pumped 
magnetometer measures the external mnaqretic field by méeking 


immersed in a nagnetic 
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Mee, OL che Sac that when an atom i 
peela, Sos erergqy levels are split. This is known as the 


Zeeman effact. Por the fields of interest, che amount of 
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the magnetic field. By measuring <2 paration between the 
zevels the magnitude of che mGq@eestc fredd can be 


det snme ned. 





This type of magnetometer usually nak 
Cesium or Rubidium vapors, or Helium gas. The cu 
ational MAD system, the AN/ASQ-81, 25 a Helium ga 
pumped magnetometer. 
In an optically pumped a2agneton 
), the sample vapor or gas (Such 2s Helium) is collected 
into an abserption cell. Gicecitasly  polarzzed Leght is 
sed through the cell giving up some of its energy in 
exciting or oumping the eléctrons of the sample gas to 
higher energy levels. These electrons then fluoresce to 
lower, metastable energy states. This is toptical pumping.' 
A detector monitors the degree of optical pumping by mneas- 


uring the transvarency of the gas cell. 
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The actual separation of the energy levels is meas- 


ured by applying a weak R. F. field which redistributes che 
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10-S nT have actuaily been measured. This SQUIDS are used 
tn both superconducting magnetometers, and superconducting 
gradiometers. 

A Josephson junction consists of a thin layer of 
insulator bétween two superconductors. This junetzon has 
meewomepesty thet a curren: can flow across it without 


developing a voltage up 


ts developed for all current 


Ga. LO eh —y algae 


In a superconducting 


mend cCOnta-ning @ peir of 


measure the 


This system is 


mm ¢ 


tions in only one file 


=O SOn@ Maxzaumn current. 


amount of magnetic 


a Vector magnetometsr 


A voltage 


Values greatsar than the maxinmuo 


mMagnezometer a superconducting 
Josephson =s used 


£1 we 


junctions 
penetrating 


which me2sures 


id component 


ato 








er 


BY uSing more than one superconducting loop, it is 
possible to measure magnetic field gradients, and construct 
a superconducting gradiometer. This is alsc a component, or 
Vector, sensor. 

These instruments have been used extensively in «he 
laboratcry and +*to some 2xtent in geophysical work, bute 
considerable engineering problems renain to be solved before 


operational Navy use can be contemolated. 


C. MAD SIGNAL AND BANDPASS 


ime source of the Signal 
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The MAD signal results fron moving a negnetometer 
through the magnetic field of a submarine, which can be 
approximated by the field of a magnetic dipole. Figure 2.4 
devicts the formation of a submarine caused anomaly. 

The magnetic moment of an object in the earth's 
magnetic field can be due +o parmanent magnetization, 
magnetization induced by the earth's eee eee O27 

C 


Gombpenation of bcth. In the f a submarine hull both 


s 5 
causes are present with a Small amount of permanent 
Magnetization produced by hull stress in metal component 
during cons*ruction and stress caused by submarine diving 
epdeeoatrcacing. This, howaver, iS @ minor contribution. The 
most importan* constituent of the submarine magnetic moment 
is due to a induced in the hull by the voresence 
of the geomagnetic field. [Ref. 21] 


The induced field of a submarine depends on che 


effective permeability along the vertical, athwartshipos, and 
longitudinal axes of the snbmarine. This infornetson , 


taken tcgether with the strength and dip anaqle of “he 
geomagnetic field, yielis a fairly precise calculation of 


the magnetic anomaly produced by che submarine. 'dDeverning' 
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of the submarine huli cancels out the permanent kull 
magnetization leaving the induced magnetization 2s the 
principal signal source. [{ Ref. 21] 

The MAD signal is approximated as the projection of 
+he submarine dipole field onto the geomagnetic field 
vector. This approximation is good because *he maaqnitucde of 
+he dipcle field is very much smaller than the maaqnitude of 
+he geomagnetic field (a few nanoTeslas as oppesed to 
approximately 50000 nT). Therefore, whenever the divole 


field is perpendicular to the earth's fleald ars 
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signal will result. The signal recorded by the AN 
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path of this sigqnal. BeicGinesi2s5 ) Guelatati 


some aspects of the MAD signal. 


2. Anderson Functions 


The submarine anomaly signal shape is a function of 
the dip angl¢ of the geomagnetic field, the magnetic 


hea 
of the eircrafit, magnetic heading of the submarine dip 


ole 
and the lateral range between the aircralt= and the subma- 
mme. Thess factors detwermine the *A’® coerficieants for the 


Andexrson functions below. 
In 1949, J. E. Anderson of the Naval Air Development 


Center determined that the MAD Signal, obtainad aiong any 


it, 


Somesce, consisted cf a lineer combination o three basic 


components. Different shaped signals could by changing <he 
Begmore2onsS: contribution of these basic components. The 


Meetete- ical revoresentation of thes? coflbonents ars r10V 
eferred to as the Anderson functions. Using the dimension- 
less parameter ‘'b' (defined as the distance traveled aloag 
mme aiccrat. track divided by the s 
Beint Of approach (CPA), b=OQ at CPA) the an 


represented by 


38 








mn re re em ee a ee 


| 


" 


ne 





Zere Signed 


| 
| 


i zere —_ 
La Sas ee giqnal "a 


per? 
—— S$: qrol 





Pos rtNe 
S tg rel Regen 


Heomagnat 
Frei e 


= Ts eS eee mS ee ee ee eee ee ee oe RR i Ee gee A ee eS eS oe Oe ee ee et ci Ee ee coi ee ce ee RD oe) ei te ei CE) ae ee 


; © Te a cm ee ee 
a 


ee ee ee ee ce ee eee =e 


Figure 2.5 Qualitative Aspects of MAD Signal 
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M A2 b©= + Al bd +A0 
Bs So a-e-- 0 25 ------------------- (2- 1a) 
By APs 
Z3 (b2 +1) 
M 
or Bs = ae (AO £0 (5) + Al Z1(b) + AZ £2¢b) ) (2-15) 
where M = ye hei ag monent of target divole 
Z = teral range between target and plane at CPA 
1 
£0(D) = seer rr er we = 
A 
{b2 + 1) 
b 
UG) tS 
a7 
(D2 + 1) 
b2 
£2(b) = wert rr = woe ---- 
a7 2 
(Dt al) 

The strength of the detected siqnal is seen to Fall 
ort as the cube of the distance between the aircraft and the 
eo get . 

By analyzing the functions and ‘their coefficients, 
i* can be shown that the optimum orientation for naxinun 
anomaly sSiaqnal datection is when the aircraft, Target 
dipole, and geomagnetic field are lined up together as 
closely as possible Sd2cliveally, this eecCurs when the 
submarine momen= and the aircraft's track are oriented 
Meetn-Souch. fRef. 14, 21] 

The Fourier transform of ‘the component functions 
were +aken to determine the frequency distribution of energy 
in the MAD signal. The ancmaly signal components are shown 
mo remiiees 2.6 through 2.8 (Anderson) [Ref. 14}. The 
PewrLTer transforms of the signal components with a platforz 
WetoCety Of 150meao=s aresshown 2n Paguce 2.9. 
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Figure 2.9 


of MAD 


Prequency Spectrum (Fourier Transform) 


Signal Components 


ran 





3. ih 
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MAD Basbter Band pass 


The MAD Signal, being a transient, is essentially a 
broadband (as opposed to a discrete freguency) signal. The 
frequency spectra of energy content shows how the signal is 
contained in a certain frequency range. In order +o screen 
out unwanted noise, bandpass filtering is used in processing 
magnetic enomaly detections. 

The major factors influencing the determination of 
an optimum MAD filter frequency ére aircraft speed and «he 
Brame cCange from the aircraft to the target at CPA. Factors 
exerting aminor influence includ2 dip angis, aircraft? 
heading, and target dipole orientation. 


The range in frequency variation is about an octave, 


with the higher values occuring when the passes are made 
parallel to the axis of the dipole and the lower when «he 
pass is perpendicular to <zhe dipoic. Anderson empirically 
determined that the center frequency of the filter is given 
by 

Oe 4 ¥ 

a a (2-7) 
Z 
where f£ = frequency in Hertz 
Vi = seivsert tere locizy 27 fi/sec 
and 2 = greatest anticipated CPA 
Tange in feet 


The) €SQ=—81 Bandpass Pilter has highpass set<irgs of 

Cee 0.06, 0.08, and 0.1 Hz, and lownvass settings of O«2, 
eieeeOsce and 2.Q0eH% (OrLOnWServicse Dagest 26) {ret. 22). 
There are no recommended sattings for normal operation since 
h 


Deackaroundmnoaseevaries, The filter characteristics for t 
Sreice =e O.omet Settings are shown in Figure 2.1 


Service Digest 28) [Ref. 23 }. 
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The figure is a representation of th2 adjustable band pass 
filter of the ASQ-81 by itself. Other parts of the ASQ-81 
system addin an additional high pass filter which adds 
another 12 dB/octave roll off to the low frequency end in 
soc Se 2.10. Thus Figure 2.10 would represent the charac- 
teristics ef the entire aAaS9-81 system if the roll off of 36 
aqdB/foctave is changed +o 48 dB/foctave. [Ref. 24] 


1.5 


1.0 


es 36 dB PER 


OCTAVE 
RELATIVE REJECTION 


AMPLITUDE 
0.1 


| 24dB PER 
OCTAVE 
| REJECTION 





0.01 


0.01 j 0.1 1.0 Vy 


FREQUENCY IN Hz 


ae 


Figure 2.10 ASQ-81 Filter Characteristics: 0.06 - 0.6 Hz 


Beang equacitOn 2-2, ame Cancer frequency tor> aj SfAD 
signal with a CPA slant range £ 300 Feet and a a 
MeLocity of 220 knots is 90.49 Bz. Using a lowpass filter 
setting above 0.6 Hz would not be very useful. 

Ore reason the ASQ-81 filter extends up co 2.0 Hz 
Meseeced Of just 0.6 Hz is the fcllowina. The sarly Service 


Test Engineering Model (STEM) ASQ-81 ‘that preceed2d «he 
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current ASQ-81 production model had two bandpass filters, 
one adjustable from .04 to .6 HZ, and the cther fixed 2% .75 
mo. 10 Hz. During tests, the .75 t59 10 Hz channel proved to 
be helpful in monitoring the STEM ASQ-81 "systen" noise 
level. During design reviews for the preduction ASQ-81, it 
was decided that it might be useful to retain the ability to 
non + or "system" Reto. wm nerOnder 25 00 thas easily, the .75 
=0 19 Hz band was dropped and a2Hz poSition was added to 


the adjustable filter. (Ref. 24] 


The center frequency at 1200 feet and 160 knots is 
moo HZ. Consequently, the highpass filter settings could 
affect slower, longer range MAD detections. 

The choice cf the highpass (lower end of bandpass) 


- 


sattings have to be made ter consijering the noise present 
2 


Q Pp 
Ith 


of systen operation 
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A. INTRODUCTION 


On a practical level MAD noise is defined as magnetic 
in the MAD ovoassband (0.04 to 2.0 


disturbances falling with 
the nommal uopec limit, 0.04 to 0.6 


Meg OL Since 0.6 Hz is 
Hz) and having an amplitude greater than 0.091 nanoTeslas 
(the system sensitivity of the AN/ASQ-31 MAD systen). 
HAD noise sources can be divided into the following 
categories: 
-~Equipment Noise 
Seeccra.t Platform Noises 


~Aircraf= Manusver Noiss 


cr 


-~Gradient Noise Due to paaea stain Hotton Seasouch the 
Gecmagnetic Fisid 

~Geologic Noise 

-~Noise om Wind Waves and Swells 


1S 2e 
~Geomagnetic Noisea 


Pee Sen SOR, PLATPORS AND MANUSVER NOTSE 


Senscr noise is the self-noise generated by the opera- 
tion of the equipment itself. Geeswican Be Beartially due <9 
She fact that the detectcr elemerc= is misaiigned with 
respect =o the geomagnetic field vecter. Changes in lamp 
intensity, vohotodetector noise and noise in the electronic 

Becutet= Can also contribute to sensor noise (Ref. 21). The 


Belt-nortse limitation of the ASO-81 is 0.01 nt. 


Placform noise is generated by components fixed to the 
oe Ca ac Z2 §the 7 Tein. 2) Or the sensor (Ref. 23). 
Permanent, induced and eddy-curzten*t magnetic fields are 
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—— a 


associated with «he airframe. Permanent magnetic fields are 


j 4 
©) 


due to aircraft structure or equipment having ferromagret 
parts. This field changes its orientation with respect to 
the geomagnetic field vector as the aircraft manuevers, 
causing field fiuctuations near the magnetomets 


PeisetoOrnm NOSSsS ts also induced ia aircraft erro 


ma 
Pemectures by the geomagnetic field. Similarly, eddy-cu 
Memcs az=e induced in aircraft skin, cribs, and frames, and 
these currents , in turn, cause additional magnetic Field 
Bails » apid aircraft manuevers will induces changes in the 
Magnetic field sensed by the magnetometer. 
Platform noise in aircraft mounted sensors is ccuntered 
by applying equal and opposite maaqnetic fields «+2 the sensor 
at 


in a process called comps: S5Onle Towed MAD systems are 


essentially free from this 


C. GRADIENT NOISE 


Gradient noise can be divided into turn noise andi neise 
due to changes in altitude 

Mien Moise is a probl2m when "WAD ‘trapping’ cr using AAD 
Eeeecrackirg a target. The earth's Magnetic Fieid has a 
horizontal gradient (in this case the magnitude varying with 


u 
Soo SKS = ; 
a» CazecCe en. Of tae 


dat iewde). bee @ hee airesatt @oves in = 
gradient the field strength chances. The noise iue +c a MAD 
€Sapping or hunting circl2 is centered in frequency at the 


Meciprocal of the tine «aken to complete one tevolution. In 
#*he case of a «=wo-ainut?2 circle, <=he noise would be centered 
oe 0.008338 Hz, well belsw the filter used in the ASQ-961. 
Set herizon-al gradient notse due to flying 2 cleverleaf 
pat clon, EOtemme emest pert also falis below the 44D 
passband. [{Ref. 21] 


Wes=-1caAil gradient noise» is du 


iD 


to chapets i: se 
moto O2:005 <<c OO.) ar 


it 


altitude. An altitude gtadient oO 





per foot exists in the earth's main field. In areas of 
geological anomalies this gradient is even larger. east 
Miewomde Of AlECraft pitch changes can c¢ ¢ a maaqnetic 


s 
be of concern. 


( 


& 
Bieia f£luctuation of sufficient amplitude to 
To avoid vertical gradient noise, altitude compensation 


equipment is used. [ Ref. 23] 


Doe GNVERONMENTAL NOISE 


Magnetic noise from sources ?2xistingq in the natural 
environment include geologic noise, temporal variazticn in 
the earth's magnetic field, and neise due *o ocean waves and 
swells. 


Geologic noise has its source in haturally occuring 
Magnetic anomalies caused bv magnetic material present in 
fame earth's crust. When the sensor passes over geological 


anomalies, the relative motion causes a MAD-like signal to 


> 
¢ 


be recorded. Geologic noise is usually more pronounced i 
Shallow water as the sensor is much closar te the source of 
the noise. Geological magnetic anomalies ar2 »sften 2ssoci- 
ated with such cceancgraphic features as seamounts and ocean 
ridges. 

Sea water is a conducting medium which is transported by 
the physical motion of water waves in *he presence 
Secmaqre ic field. This motion induces currents in 


These currents give rise ts sscondary magnetic fields, which 


eea@ weetOrially to the quasistatic, qeomagnetic £1211 


= 


(Weaver) [Ref. 25]. These fields can be detected at siqni 


tty 
j 


icant distances above the sea surface and fall off exponena- 
moe lly @eth altitwees Pajew@ee 36125 a8 plot ror several 
Surface wave periods of the induced magnetic field oer nete 

amplitude of the surface wave. These induced fields can be 


a problem at the low altitudes where 4AD is used. 
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aye 








Ee GEOMAGNETIC NOISE 


Temporal variations in *he earth's magnetic field with 
frequencies in the MAD bandpass and amplitudes greater than 
0.01 nT have become known aS geomagnetic nois= in +*the MAD 
literature. 

Quiet daily variations, such as the Sq and L v 
have periods sufficientiy long +9 fall far below the MAD 


passband. 


Geomaqnetic storms have been discussed previously. 
Rapid fluctuations with high amplitude and falling within 
the MAD passband occur in connection with geomagnetic 


StOIms. 
Geomaqnetic microrpulsation Conon te the last category 


of geomagretic noise to be discussed. 


Geomagnetic microoulsations are tapid fluctuations 
of the earth's magnetic field with periods from 0.2 seconds 
to 10 minutes and amplitudes from about 0.1 nT to as high as 
aeew Tons Or NTs. These fluctuations are caused by elsc- 
+romacqnetic perturbations propagating in the magnetosphers 
as hydromagnetic waves (Nishida). [Ref. 26] 

Micropulsations are classifis py morphology, hat 


Hop OY =xXaMining perieds, amplitudes, times of cccurence and 


other observed characteristics. Mtcropulsations have been 
placed into two broad categories: irreqular and continuous. 
The irreqular pulsaticns are represented by the symbol Pi, 
and the continucus pulsations by Pc. Table II shows the 


breakdown by period of the various nicropulsation classes. 


iD 


It can be seen that for MAD geonagnetic noise th 
fe) 


th 


§ 


t 


MeowSe “frequency Poel, Pc2, Pc3, and P21 pulsations are 


Seneerest. 
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YABLE If 


Geomagnetic Micropuisation Classes 


NOTATION PERIOD (sec) AVERAGE AMPLITUDE (nT) 

Pc} Ore2 = 5.0 9.05 a | 

Pc2 >” 19 On 1 =. 0 

Pc3 10 - 45 0.1 - 1.90 

Pcy G5 - 150 Os.1 =~ 10 

Pc) oe sb O0 ia) 410 

ea 1 - 4Q OOS ela 

Pi2 GO = 159 dugee 

(arteter Jacobs) 


Pol pulsations are regular sinusoidal oscillations 


e 
With periods normalilly falling in the 0.3 to 4 second (0.25 
memos sos fz cCrequency) rang 1 

bursts and gradually develop into 2 series of puisations 
which could last for hours. They may also occur as consécu- 
tive groups of pulsations with sharply varying freguencies. 
The average amplitude of Pci's 1s 9.05 to 0.1 nT and they 
jena tO have a single well defined frequency. Pct vulsa- 
tions with frequency less than 0.5 HZ aze more commen ar 
Monee sea titudes than at mid latitudes. Thess oulsacions 
Oeecur in the daytime in the auroral zone, and at nigh« and 
@arly morning hours in lower latituies. Pcits are charac- 
tes of the 


Becastic Of the quiet and weakly disturbed st 


é 
geomagnetic field and show an incre Ley wOn= =o 
nd 


a 
two hours oefore, and tour to seven days after, magnetic 
eocm. {Ret. 6 } 


Pez,eavd. hes Dilsa=2 ons ace rouped t 
2 


fv 


QO 
@earacteristics. The amplitudes of eectzulet.ons are 
hg 


O 
usuelly under 9.5 nT and the *«ypoical frequenc ange is 0.03 
mee 0.2 HZ. These pulsations are normally found during the 
meee ten ciel = a2ce2v¥2 cy Foacnutng a Maximum arcund noon. ?Pc3 
piesactltOhe Sne@wW 2 seasonal variation with a minigun of 


SeesvVity occulriag during winter. 


a 





Pc4 and Pc5 are large amplitude fluctuations, but 
fall below the frequency band of interest for MAD. 

Pi1 pulsations have an irregular form with an 
average amplitude of 0.071 to 0.1 nT and a freguency mainly 
fmoecrne O10 to 0.17 Hz range. Spectral analysis of these 
pulsations show a braod band of frequencies. Pil amplitudes 
have maximum values in the auroral zones, with the intensity 
of the puisations decreasing with decreasing latitude. 
Piil's are normally cbserved in the late night and early 
morning hours, and show an incr3aase in activity with 
increased geomagnetic field disturbance. [Ref. 6] 

G2omagnetic micropulsations can be observed anywhers 
on the globe, at various times of day and year, and in both 


quiet and disturbed geomagnetic field conditions. 
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A. INTRODUCTION 


Previous chapters have defined ae tte activity as 
s chapt 


memecppltes =O Magnetic Anomaly Detecticn. 32 te lela 


e ia 
methods of e¢valuating that “oe Willi Ge exanised, 


including methods currently in use by the fleet. 


Eee GEOMAGNETIC INDICES 


A gecmaqnetic index is Simply a ne 
oyna 


ard describe time variations 


Becuiting from solar-terr 


rc 
indices are commonly used +O eExpre 
a 


depict the character of gacmagnetic 


day. 


mor  =he most part, ces developed as 


Tenge indices, measuring + ove fhweh and 


ct 
th 
{~ 
‘D 
> Fe 
fu 
a 


low values for differen: measured during 


the day by magnetic observatories (Lincoln) [Ref. 27]. Most 


G@erent indices are of the <cang2 types, but other indices 
have been developed which arte more subjective or qualitative 
wm nature. 

Geomagnetic indices are designated by a letter code such 
as: Gye, Cpr eo, CO, at Nneastes, Det, KP Ks, Kp, ak, Ak, 


ap, and Ap. There are additonal indices in use. 


ov) 


ThesGy Gi, Gp, and C9 indaees are daily magnetic fiel 


° 
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character figures. Tics Caimdexaeis tieeadaily 


Mm Oa 
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pu 
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fu ju & 
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figure for a singie observatory. T2, *h2sseal> 


cates a quiet dav, C=1a noderately iisturbed day, 
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heavily distuzbed» day. The daily international 














(D 


meqgure,)o., 25 the arithmetic pean of the C indices reported 
by participating observatories arecuni ‘the world. Gage “=hs 
deily planetary character figure, is similar to Ci, ¢xcep 
that it is derived from the values of Kp and ap. c9 is 
@emecacted Scale for Ci ancCp with single digit values 
Funning from 0 to 9 (Bartels). [Ref. 28]. 

The Q and R indices 2re quarter-hourly and hourly range 
indices respectively, taken at high latitude stations only. 

The W measure is an index of the equatorial electrojet. 
Dst is a measure of Ting current e2ffect. They are both 
amplitude indices. [Ref. 27] 


Als K, 


ine keno, AD, 2%, sake, aby and Ap indices comprise a 
group of related 33-hour range indices. Tre A. dedex is a 
Singie station cede using a gquasi-logarithmic scaie ‘from 0 
to 9 *o measure geomagnetic activity. The value of K is 


determined by first determining the difference? between «the 


lowest and highest deviations from the regular daily varia- 
tion (Sq) during a 33-hour period. Uns ce ses. (in Tee is 
GGnvezted *c the K scale based on the histcrical activity 
anges at the particular observatorv involved [Ref. 27]. 
The conversion for the Fredericksburg, Virginia observatcry 
w= given in-fTabbe III. Thal conWers@on cam aiso be applied 


men eme USAY/NOAM Observatory in Boulder, Colorado fRef. 29]. 


The KS 2ndex 25S a Standardized K iniex which is 


th 


moa ->Omn © local vartaztons and is then used +o de*et&ins 
fe) planetary 3-hour index, Kp. 


THE equivalent three=hour-range, ak, is a convetsion 


ect the K index as shown in Table IV. En™order vor detente 
mie uneicts wot ak for 2 vartiicular observatory, d@vaede the 
Meyer cange lamat of K=9 by 259. Thws tor rredericksburg 


and S3outider, ak is in 2-<nI units. 


on 
oO) 








Pesos LT 


Conversion from Range to K for Predericksburg, Va. 


K Range(nT) 

0) Oo = é 

1 5 AD 

2 10 - 20 

3 20 =- 4Q 

” 49 ~- 70 

5 Oe || 218 

6 ioe ==200 

7 2900 =- 3390 

8 S30) Sse ae 

9 > 500 

fette> penicolin) 

Pasko. Lv 


Equivalent Range ek for Given K 


3 U 5 6 7 8 S, 
ak OQ 3 7 15 27 48 80 140 240 400 
(ef-onTirncoln) 
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is the equivalent daily amplitude and 
= 


= 
average of the eight daily ak values at e@ particular observ- 
O 


17) 


peo sy . Sh Jala index is promulgated using «he name 
observatory, the Ak index fcr Fredericksburg is known as <n 
Pete tredericksburg cr A-Fred index. 

The equivalent planetary amplitude, ap, Es deter- 


eee tecommthe Kom index in a fashion similar te ‘hat of 
Ln 


ja 
(D 
tA 
ih 
O 
ie 


ct 


x 
determining ak from the K indices. The aight ap val 
@given day can then be averaged into the daily equival 
planetary amplitude Ap. These two indices are given in 2-aT 


Una tS. 





Sa GLO AGHETIC@MDICES IN BLEEBT USS 


Geomagnetic indices of interest in connection with 
MAD Operations are the K, Ak, and Ap indices. 

Fleet operators utilize tke Neen da eee oO Sk EOL 
perdicting geomagnetic activity over the antizre world 
eet. 9S0 j. This index is promulagated by the Flsaet 
Numerical Ocenaographic Center, Montarev, California, in «the 


environmental briefings received by aircraw personnel. This 


index is the Ap index as sent out from the Svace 
Environmental Services Center, Boulder, Colerado in tha 
Joint USAF/NOAA Report of Solar and Geophysical Activity 
Ret. 31). 


Tre BOUlIGer K index is available +t5 interested 
parties by telephone recording ard in the WWV and WNHVH ca 
meeeacests and is therefore Avesta obe wre O0 thoeet, Users 


fRef. 32). 


The Ak index for the Frsdericksburg, Witla oa: 1a 
ebservatory has been used in studies of gecmagnetic activity 


aS applied tc Magnetic Anomaly Detection [ Ref. 33]. 


Oe encomecical Applicability of # and K indices to MAD 


With the widest useful filter settings, «he MAD 
Demepass tanges from 0.04 to 0.6 4z (1.7 to 25 seconds ia 
period). is SCH, cneObdenr for a geomaqnetic index to he 


directly applicable for MAD use, heeswould we sensitive tc 


4 


that frequency range. 

The K indices, and the K-derived A indices, are not 
especially sensitive to the MAD range. Mavaud [ Ref. 34] 
indicates that these indices are mainly sensizive to fluctu- 
ations whose periods ar2 much lo 
the frequency range analyz2d, that is, a frequency corre- 

es 


Sponding to a period of 45 minut 
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OmiemeGeasOn form. <ass lack of sensitivity for MAD 
bandpass geomagnétic nois® is that the amplitude of gecmag- 
netic fluctuations varies inversely with frequency, so that 
the amplitude of the fluctuation increases as the as che 
frequency decreases. It can therefore be seen that the 
MiectMations with periods of an hour or great 
determine the variation range used t5 calcuiat+e the K index. 
Mae activity driving the K and A indices, will, because of 
band pass filtering, not even be observed by ch 
and the activity of interest to MAD night not 
Sere A indices at all. It can be concluded that * 
direct physical iink between MAD geomagnetic noise and 
either the K or A indices. [Ref. 35] 


3. Experimental Correlation of A and K indices with MAD 


ae ASQ~-10A Study 


Brennan and Smits [Ref. 33] fcund that gecmag- 
netic micropulsation activity was reccerded at their ASO-10A 
MAD magnetometer site in Marviand, when the A-Fredericksburg 
index was greater than 25. This cccenred everytime they were 
Meee-dang Gata With A-Fred greater than 25. It is important 
Wemgoee «iat <=hey observed additional activi~y during some 
periods when A-Fred was less than 25. 

The Brenran and Smits study tenis to vaiida 
the use of the A indices as at least qualitative indicetions 
of geomagnetic noise in the MAD cated Som: Theiz study, 
however, waS specific t> tne ASO-10A wmeaegneton 
which has a sensitivity of 0.1 nf as opposed ts the 0.01 aT 
sensitivity of the ASQ-81 systen. Time eefect of the ASQ-19 
Sensitivity is to filter out most Pct pulsations. As 2c1 
pulsations do not correlate well with «the A and K indices, 


these voulsations would «2nd to increase 


ae 


th 


mn= t2beerzag out o 








the reliability of the A-Fred or Ap indices as measures of 
MAD geomagnetic noise. Mason [Ref. 36] stated that the 
toccurence of Pcl is well known to be associated with low Kp 
values." It should be noted that an operational drawback of 
Mmimtaicering Out Gf Pel pulsations by the ASQ-10A is “hax 
the system also filters out valid signals of less than 0.1 
me amplitude. 

While information has been presented that 
Suggests the the A-Fred index can be useful for MAD gecmag- 
netic noise evaluation for the AN/ASJ-10A Systen, sufficient 
data was not presented to draw conclusions for index usage 
with the AN/ASQ-81 systen. 


Deano S| seudy 


For two weeks in April 1976, Naval Aiz 
Development Center personnel operated a geomagnetic observa- 
tory at the Atlantic Undersea Test and Evaluation Center in 
the Bahamas. The primary magnetometer used for this observ- 
atory was the ASQ-81 nagnstometer. One purpose of this 
observatory was to conmpar2 the K ind3x (as cdetermined by the 
San Juan Observatory) with geomagnetic activicy in «he MAD 
band pass. The conelusion of this stud¥ Was that the K-San 
Juan index did not correlate with geomagnetic noise in “he 
MED band pass. {eRef. 37 j 

We have undertaken a correlation analysis using 
Se fasner Z transtepmation [R@ef. 38] of *h= WADC daca 

; A 


misdished tO us by OcGhadlick. The data used is list 
Tables V and VI. 
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being from 0.14 to to 0.52. Sample size was 86. This indi- 
cates that there is at most a weak correlation between the 
observed data and the K-San Juan index. A much greater 
correlation would be raquired for the K index to be of any 
significant value for MAD operational use. 

The coefficient of correlation has values fron 
oe tO.|CU1 0. Reeve ie 55 =-1.0 or .O0@indicates perfeet 
negative or positive correlation, respectively. A valne of 
zero Signifies no correlation at all. 

A similar correlation analysis against the Ap 
index was conducted. The sample correlaticn coasfficient was 
0.51 with the .95 confidence interval for the actual coeftfi- 
eee Oc COLrelation being from 0.32 to 0.65. Again, <has 
Signifies that only a weak correlation was observed. 

The weak correlation between observed gecmag- 


netic activity and the K-San Juan index suggé 


{Nn 
cf 
in 
cf 
oa Hee J 
$v 
c+ 
cf 
; 
t- 
W 


index wceuld not be very useful in describing MAD geo 
Momse, aS this K index reflectad activity similar to that in 
the band of interes* only about one-third of ths tin 

Beomgem Correlation of the Ap index irdicates that it 
reflected activity similar to MAD ysomagneti noise 4 
one-half of the time. This is still not a very good indica- 


peeemee wid 15 Going on in the MAD band vass. 


Cc. Power Spectral Density Evaluation 


Wseeeasc of OMdoing rescarch at the Naval 
Postgraduat? School, geomagnetic activity data in the range 
Oz cre SAD band pass has been coliected and enalyzed in <=he 


form of power spectral density (PSD) curves. 

A measure or MAD band activity 
oped by integrating under the PSD curve and using compre 
mana bess Lligitts of 90.05 and 1.0 Hz. Theis a 


discussed in depth later. 
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Two sets of data ware analyzed. One set was 
taken using an induction coil to measure the fluctuations in 
one direction during May, 1980 [Ref. 39] and August to 
October, 1981 [Ref. 40]. The other data sat was collected 
uSing a Cesium vapor toatal field nagnetometer during the 
July to October, 1980 period [{Ref. 41}. 

Correlation analysis of the single-coil system 
data yielded a sample correlation coefficient of 0.401 for 
the K-Fredericksburg index, 0.780 for the A-Predericksburg 
eigex and 0.155 for the Ap index. Sample size was 9. The 
memGeieraence Gntervels for the correlation cosfficient for 
eee e, F-Fred, and Ap were -.36 to 2.85, -.55 to .76, ani 


oan =O .75, respectively. treetaea Goer thas test is 
presented in Table VII. 
EASE e Vil 
Sanql.e  Golleres Noise Data end Indices 
RMS Noise 
Date tS) Amol itude (nT) RoeeeG GA rae ap 
1 May 80 180 O20 3 10 
al aud 81 2045 GE0S 2 lB qu 
Z> hug 31 0235 Opa 3 tbs. 16 
25 Aug 81 0745 0.19 3 15 16 
Za Aua 81 1318 0.05 3 15 16 
4 Oct 31 1830 0.06 3 i 17 
Sagc- 31 0785 Ore 2 1 7 2 
Smee: 81 0630 0.03 5 7 7 
peoees, & 1 1530 Or. 02 6) i 7 
The Cesium vapor magnetometer data yielded 
Beemepbre ccerric:Snts of correlation of 0.552, (2574, “enc 
0.444 for the K-Fred, A-Fred, and Ap indices respectively. 
The sample size was 14. The .95 confiderce inte ls were 
Bee Rh -Frred from .03 £0 .d8%, A-Fred from ~.20 to .76, and for 
Pees emee —-- 12°20 .79- The Cs vapor aqnétometer data is 


orssented is Table VIII. 








aco ey I 7 


GS) Vapor RMS Noise Data and Indices, Jul-Oct, 1980 


MS Noise 

Date ane ue io ease (nt) K-Fred A-~Fred He 
25 Jul ASO 133 0 9. 10 6 onl 4 
20 Jul 1450-1520 0.06 @ 7 8 
30 Jul 1520-1551 Ore.0n. 2 Fr, 8 
al Jul 14 50-1520 0205 1 6 6 
31 Jul 1520-1555 0.04 1 6 6 
16 Aug 0030-0106 0.05 2 19 26 
16 Aug 0900-0930 0.04 2 19 26 
18 Oct 02 30-9300 0.06 3 17 ida, 
3 Oct 0300-0330 0.06 4 ly 21 
18 Oct 9630-0709 0.04 Q 17 21 
ie Oct 1030-1100 0. 04 3 a? 21 
mo Oct 1100-11390 ®.05 3 7 24 
io Oct 12 30-1300 0.04 3 7 2 
fe Oct 1290-1330 0.04 3 V7 21 

Although che sample sizes used were 00 smail to 


Stee any Meaningrul conclusions, ther2 is littl? ev 
Suggest that any of the K-Fred, A-Fred, or Ap indices 
very accurate measure of geomagnetic noise in =he ASQ-31 AAD 


band pass. 


me CGomeclatszon Conclusions 


Although some weak correlation does exist 
between the K indices, <*+hs A-Fred index, the Ap index and 


geomagnetic noise in «he MA 

Micedental and ind@isect, dein s 
Der weer ths activity in the MAD band and in the iower 
Be=GQUency activity that influences the &K and & indices. 
Mieme 27pe2ces acs not dirsctly influenced bv aczivity in the 
MAD pand. The correlation that does sxist does not appsar +o 
be sufficiently high to enable these indices *9 yieid accu- 
rate indications of the actual MAD band activity. The use 
Oz these indices for anything except the reughes*t quaiita- 
sa¥e@ @Stimation of achive in the MAD band pass is a20*% 


recommended. 


oo 








D. PROPOSED GEOMAGNETIC INDICES FOR MAD 


Overail geomagnetic acitivity is analyzed in both «he 
time and frequency domains. Geomagnetic noise indices for 
the MAD band pass could be developed in either of these 
domains. 


The spatial coheren zi MAD geomagnetic noise has not 


ne C 

yet been adequately determined. This information would be 

necessary in order to determine the number and locaticn of 
Pp 


mini-observatories for an sperational MAD noise index. 


One way to develop an index of MAD g20magnetic 
activity would be to establish nini-observatories near bases 


rom which MAD operations areée conducted. These observato- 


tt 


ries would use ASO-81 magnetometers or different magnstome- 
ters with ASQ-81 filter networks, and could in real time 
record the geomagnetic noise in the MAD Ddand. A neasure 
such as the maximum peak-to-peak (9r possibly the average 
peak-to-peak) noise ina civen time period couid then be 
disseminated to flight crews operating in the area covered 
by that index. Obviously, the spati21l ccherence of MAD band 
Seravity is important in making such a2 system work. bh oe 
+ype of mini-observatory has been suggested by References 31 
aga 35. 


2- Frequency Domain Inde 


Present fleet procedures examine MAD noise such as 


3 


System and manuever nois®? in terms of the amplitude of the 
2 


Piien= ve liom {| Ret. 30]. an index of geoMagnatic noise in the 


MAD band pass wouid therefore bea of greatest usezulness to 
teweec. Operator LF it Were in units of the amplitude of 
she signal as seen by the MAD eguipmenct. 








The method proposed to derive a AM 
the frequency domain begins with obt h 
tral density of the activity in th2 MAD band by Fo 
analysis Cf ‘the tine series jata D 
magnetometer. 

By intergrating under the PSD curve over the limits 


of the MAD bandpass, a value inunits of amplitude? will 
result. Taking the squar2 roet cf this value will yield an 
RMS amplitude. Equation 4-1 represents the derivation of 


this index. 


u ree 
MAD Index = if PSD(E) dt (4-1) 
where MAD Index is in nT (gammas) 
£f = frequency (Hz) "Tk 
u = upper bandpass limit 
i = lower bandpass limit 
and PSD(£) = power spectral density 
(nP27 nz) 
Ge (Ghatecteriscics Of the HAD Elster Eom er 


not being an ideal bandpass fiiter) could be applied prioz 
momene in seqretion. The integration itself could be done by 
Peenef ad peoznat—by-volint numerical integration or by first 


modelling the PSD curve by polynomial curve fitting and then 


integrating the polynomial over the range of the band pass. 
iota anticipated thet this couid b2 done in ciose to real 
time by a digital computer, possibly by a desk tcp computer 


such as the HP9845. 

Pion yoe Ore S2nseeeuutzi2zei could be the ASQ-81 
magrnetcmeter, other total field magnetometer, or pessibly an 
ortthogenal 3-coii system whose signals can be combined to 
veela hey projection on the total Field vwecter of «he fluc- 


Beeceson OF 


j-- 
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tWaticns. A Singl2= coil system oriented 


D 


used. 


1D 


th 
the Gatrrzh's magnetic field vector could also b 











A three coil system which is used t9 yield an RMS 
amplitude is currently in research use at the Naval 


Postgraduate School. 


3. Predictions cf Geomagnetic Activity 


The proposed indices discussed above are intended to 
be real time measures of tne aqeomagnetic noise in =he MAD 
bandpass. Whether or not such activity can be predicted 
ahead of time needs to be looked into. 

While there is no model for the background component 
of geomagnetic noise, work has been done on estimating the 
PuGure actaivaty Cf micropulsations, nctably by Fraser-Smit 
Sn the case of Pct pulsations [Ref. 42, 43}. By extendin 
She prediction technique for Pcl pulsations to the Pc2, Pc3, 
and Pil pulsations, the occurrence of geomagnetic micropul- 
Seeeo2em activity in the MAD band might b= predicted. 
Gembpanang this prediction with sreal-time solar flare infor- 
Mation should give the capability +9 disseminate real-*ime 


and estimated futute MAD index values to fleet users. 
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Previous studies have le nor rT conclusion that 


d 
currently used gzomagnetic indices are not accurate measu 


res 
of geomagnetic noise in the MAD band pass. Experimental 
equipment has been utilized, and computer software written 
mammmomd== +o confirm <tihis conclusion, and to develop a2 


replacemen> means of evaluating MAD geomagnetic noise. 


A. EQUIPMENT CONF IGUSATION 


Experimental equipment, acquired as varz of the Naval 
Pestgqraduat2 School gqeomagnestics research proaqram, has been 
utilized in *he 2ffort to develop a usabls MAD index. Ths 
Peete 25 If WSe “hn Other projects of the gecnagnetics 

sociated equipmen= are set 
2 


research greup. The sensors and as 
) 


memtert remete site operati mimcnitoring and data 
analysis located at tha Naval PDostraraduaces SG HOO 
Mesertpu20ns Or the data collection system and data aralysis 


The data acquisition system illustrated in 


rr 
H é 
\Q 
j= 
4 
iD 


5.1 reveals the following majcr components: 
~cceil antenna sensors (3) 
Be caipwie.ems (3) 
-~Signal conditioners {amplifiers) (3) 
-pulse code modulation system (1) 
~radio «ransmitter (1) 
-radioc receiver (1) 


~instzrumentaticn tape recerder (1) 


ee, 
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2nd Coil Antenne Sensers 


Three coil sensors ere used in this 
Bach sensor is a self-supporting, continuously woun 
center-tapoed coil antennas manufactured by Elma Eng 
Peto Al+o, Californie, Peom about 5460 turns oo 
copper magnet wire. The colls weigh approximately 50 kg ¢ach 
with dimensions as depicted in Figure 5.2. The dimensions 


or the sensor are constrained by e dimensions of the 


}~ 


el 
largest glass sphere that 1s commercial 
< aL 


20 Shams 2nd i2s self- 


y available. These 
S 


Spheres are used to enclose the 


Q) 


Guring und2rweter 


O° 
experiments. The coil résistance is 1 


inductance is approximately 9.31 henries. The < 


et 
ara mounted orthogonaly on a nonmagnetic frame (Fisure 5.3 
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Pigure 5.2 Sensor Dinensions 


The preamplifier used w2s the model 13-10A low 
Dp n 
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ee weren=eqd Dyept. Alan Paitllizros of SRE 
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Figure 5.3 Sensor Mounting Block 


jMmoeeamemoneal. The ftEtanal stage of the amplifier contains an 
feeeawe Jowepass Sle er which orowides a sharp cutofé for 
y= 
a 


frequencies 2bove 20 Hz. The cya 
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MapiersmCe bese cada 2.65 MiLiivolts is 60 «dB. 
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Gc. Sigma beComditio ners 


The signal conditioners receive the analog 
Signals from ‘the coil preamplifiers, amplify them on «he 


order of 30 4B, and limits signals with peak amplitudes of 


7.5 volts from entering the pulse coi2e modulation sys*en. 
d. Pulse Code Modulation (PCM) Systen 


The pulse code modulation system chosen for use 
is one designed and manufactured by Or. Robert Loewe, 
bemiecom, Inc. The system features 15 channel anaiog input 
capability and offers selectable sampling rages of 2, 4, 8, 
16, 32, 64, and 128 samples per second. By appropriately 
jumpering the analog input pins, the samplin rate mav be 
meemeesed by a factcr of 5. 

The PCM syst? 


lator, and associated Cif 
a 


len meomnotidess) 25Ceystal oscil- 

ome heeGeatea Gc. =Ctrtary tc develop 
the clocking puises, TiecmmOmenetn-meenOSs analog muiti- 
plexer, a 16 channel, fo Dit netos andaregs £O digital 
converter and associated circuitry *o provide the pulse 
Coding. ieee cael check  Oscitliazor operating at a 
zai 


alo 


x 
ip 


Frequency of 24.576 kHz produces a sauare wave output 


loss rate of 7 bit in 106. The clock pulses gate the 


\Q 


p 
rs 


fu 


Multiplexer, aneweq “csumdagd22al COnwWerte> and associate 


col log cn Secuctry chat LOrm “the pulse code words. By 


(D 


basic output is a Bi-phased pulse coded siaral. 

The data is organized in frames. Each frame is 
headed by a sync code word which is followed sequ 
the pulse coded samples from PCM channels 1 throu 
syne code word is a pulse coded digital word wit 
value between QO and 4095. Mies Vosa 2S Dl5es 
Per dwiamedeon theecitcuat boards Thais code werd is essential 


to the decoding process. 
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im <<ne swat ze Operation of this systen, a 
sampling rate of 32 samples per second was utilized. Only 
one sample per coil per frame was analyzed, making use of 
PCM chanrels 2, 3, and 4 only. 


éeé. Transmission and Recording 


After the data has beer PCM encoded, it is 
+ransmitted by a VHF radio link back +o 4 receiver Located 
at the Naval Postgraduate School, where it is currently 
recorded on an instrumentation tape record=c fer iéter 


analysis. 


ct 


2. Data Analysis Eguipmen 





Currently, the recorded PCM data is played back into 
a PCM decoder and associated equipment which g 


€ a 
nine-track 800 BPI computer tape containing the decoded 


sensor Gata. Thas computer tape is then input intc the 
Naval Postgraduate School IBM 3033 compoutez for analysis. It 


*s in the computer software that the sensor s 
Zunctior is applied, spectral anaivsis perfor 
MAD index generated. 


Ee DADA ANALYSIS SOFTWARE 


Hy 


AS was noted ecariier, a mainframe computer was utilized 


to perrozm the spectral an peEY =ne Fz ansre> Ewunc- 
<t20n, convert to power spe y 


Cc 1 
generate tne RMS noise anolitude MA 


Epa (Using equaticr 
Crt) . Pics PE@etam is OWretten in FORTRAN IV and is 
discussed in brief below. Mi copy of the oregram canbe 


found in Appendix A. 
The, Walter peecram is divided ints sections which verforn 
Beef eCliowing funct=ons: 
~Data input 


-~Fourier analysis o= tim2= series data 











-Application of system transfer function 

-Projection of field components onto total field vecter 
-Data averaging 

-Curve fitting and calculation of MAD index 


-Calculation and plotting of power svéctral density 
fe Data Input 


Data input is accsmplished with the aid of 4 subrou- 
tine package supplied by “eee Scan-On cof the Naval 
Postgraduate School Department of Oceancgraphy. His subrou- 
*#2ne (Calied "SUBROUTINE RD‘) serves as a FORTEAN 'READ!® 
statement, takin the PCM data off the computer tape and 
converting it into integer format with a value between 0 and 
#095. The data input section of the main program takes this 
M@mceg@ger value and converts it to a "REAL' number ard normali- 
izes it to represent a voltage value between <-5.0 and +5.9 
volts. This section also sorts the input data matching the 
PCM channel to the data array representing “ths appropriate 


@oul. 


2, Fourier Anaiy 


Ie. 
In 


Thies @ig]e series data 15 next cenverted +o he 
EtTequency domain by utilizing a subroutine (called 'FOURT') 
which performs a Fast Fourier Iransforn (FFT). The subrou- 
tine is one available <9 users at the Waval Postqraduate 
Bemool and weilizes the Coolsy-Tukey FFT algorithm. Yrurther 
information about this subroutine can be fronnd in “the 


Paogram iisting in Appendix A. 


Pee eri iGeraon Cf Tmanst@®e Function and Total Field 
DE O1eC teen 


Tre neke s@ction of code wapplies the system transter 


PimectuZons LOrethnewthree coils ©to tha frequency domain data. 


das 





The transfer functions are given as straight iine segments 
which were found by least-squares approximation. The data 
enters this section in amplitude units of volts, ard is 


converted into nanoTesias by the transfer function. 
Following the application of the system trarnsfe 
mameticn to the coil data, the program next calculate h 
Bojecticn of this data onto the earth's magnetic field 
Weetor (total field projection). Tags is done by fers 
applying the local ae? variation (declination) to the 
North-South (X) an 


Mine the horizont 


d East-West (Y) coll information to deter- 
Pe tieclda Gcoupeomen:. “The lecal magnetic dio 


a 
(Qnelination angl is then used to project the vertica 


N 


at 
[—co2l) and korizontal field fluctuations conto the total 
On, 


mela f£luctuatio 
4. Data Averaging 


The previous program sections exist insid2 of 2 
do-loop which enables the anaiysis of a long peric 
without a prohibitive need for storaqe space. T 


fe) 
includes accumulator arrays for each field component ana the 


total field projection. Ipc sec wation data is converzed 
[PaeopOWer prior to storage. Tiuswes.s lone by taking =he 
Magnatuce of the fluctuztion component, adayideng by) she 
number of sample points, and then squaring the vaiue. After 
the pregqram passes through *the averaging loop icr the last 
Gomi, tae arithmetic average is ‘taken for each frequency 
Peent Cn the arrays. At his stag? the power spectrum is 
Multiplied by the sample period to datermine power spectral 
density. 
5- MAD Index Caiculation 


The next section cecmputes the RMS MAD rnois? index 


previousiy discussed. A polynomial curve fi%* is performed on 
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the total field PSD using an available subroutine ('CHBFT"'). 
The resulting polynomial is then integrated over the limits 


of the various ASQ-81 band pass settings. 
Oe PlOtrrapomot Pomsr Spectral Density 


Plots of the power spectral density of each of thse 
field components and the total fie1l] projection are gener- 
ated in the last section of the progran. Pigs) 2s dene by 
converting the fluctuation power spectral density tec deci- 
bels (dB) referenced +o 1 nanoTesla® per Hertz. A Véersatec 
Peotcting subroutin ("PLOITP') ais then calied +o actually 


generate the plots. 


Sent TIAL SYSTEM OPERATION 


The NPS MAD index system was initially placed into cper- 
Memoteween oh] Coil sensors located in the La Mesa Village 
housing area near the Naval Postgraduate School, MWYonterey, 
California. System checkout was accomplished in June, 1982. 
Gmeewtull system wes placed into opsreation on 25 July 1382 
and 18 August 1982 in conjunction with similar measurements 
taken or the floor cf Monterey Bay. 

hmiemwMaAp andex output and power soecteal density plots of 
Siemecrtel fzeld FEluctuatcion for 25 July 1982, 1237-1406 
Weeceal (2037-22062), and Wee AucGte.e 1962, 91271-0250 local 
(0921-10502) and 0507-0635 local (1307-14362) are shown in 
Tables MIX, Rie Xen 19d Figukss 5.4, SA SOO ace 


respectively. 


ee 





TASLE IX 


en 


Monterey 7 


Zein 2a 


MAD Noise 2037-22062, 


NanoT2sla 


1 


1 Ganma 


0.19083 GAMMAS 
O205 ero 0.22) HERTZ 


ws NS 
Sa a 
0G ‘aa 
WRI WMI 
eo a oi 
¥; Fa 
=a) > 0k a) 
mMCt To 
(D e« We 
© oS 
GO oy 
as WO 
=O eG 
eet oe Et 
e a 
Os Ost 
an) Leis 
e e 
Sa | © 
It It 
4 ee td @ 
hwM mM 
aAW AM 
A= ar 
HA, HAs 
Cy A 
Ae AR 
Qo «Cert 
sm 2:0 


; 


HERTZ 


0.0991 GAMMAS 
9.69 


0.06 TO 


tI Ns 


MAS 
817558 
A 
HER 


Alt 
4 Q 
AM 
60 


0.0708 G 

Oa (0R 1G 70. 
Q 
at 


MAD LNDEX= 
BANDPASS: 
De 
De ASS 





CO 
se 





0°S Sh O°h G°e ‘One G°Z ie 
(ZH) AONINOAUI 


IWOOT S9OHT-LEZT SZ286T INS SZ 
SOVITIA VSAN VI SNOILOGLONd AIdId 
WLOL JO ALISNIG IWYLOGIdS NIMOd 


S° 


Oh- 


02 


0¢ 


Ot 


09 


08 


Gi) Sieaqlora 


Wee Ts 
H/,LlU 1 3 


(2 


La Mesa Village 


PSP wD) JUL 82, 2037-22062, 


Pigure 5.4 


1S 





i; 





TABLE X 


tareay, CA 


fon 


18 AUG 82, 


MAD Noise 0921-10502, 


1 NanoTesla 


Gamna 


0.3644 GAMHAS 
OnctmetOn 0az0 Fant Z 


MAD INDEX= 
BAN UP ASS : 


80 


ee a ee 





0°s S*h O'h Sue a. € S'% 
(ZH) AONANDAYS 


IWOOT O0SZ0-TZTO ‘Z86T ONV BI 
JOVITIA VSIW VI SNOILOGLONd AIGI4 
IVLOL JO ALISNIG IVYLOIdS YIMOd 


GS, U 


mo. € 


Sa0 


09> 


Oi 


Qc- 


=) 
Zaye Eee) Siddtosd 


© 
NS 


Oh 


09 


La Mesa Village 


PS® 18 BUG 62, 0921-10502, 


Te) 
® 
wn 
a} 
‘4 
ba 
Lo A) 
“rf 
fe, 


81 








TABLE XT 


CA 


Nome okey, 


18 AUG 32, 


MAD Noise 1307-14362, 


1 NanoTesla 


1 Gamma 


GAMNMAS 
2200 PERIZ 


0. 3504 


0.04 TO 


82 








Sth O'h S°E€ O°€ 
(ZH) XAONTNODYA 


IWOOT 9€90-L0S0 ‘Z86T ONV BT 
FTOVITIA VSAN VI SNOILOGLONd AISII 
IWLOL JO ALISNIG IVYLIIdS YIMOd 


mie 


S& 
™N 


(ZH/,1U T 424) STaaIoad 
1307-14362, La Mesa Village 


83 


PSD 18 AUG 82, 


Pagquse 5.6 











em CONCLUSION 


Aralysis of available information regarding current 
geomagnetic indices and the actual level cf geomagnetic 
noise in the MAD bandpass indicates that the currently used 
indices, the K and Ap indices, are not valid for MAD opera- 
on S. It is therefore desirable to derive a new index 
which more accurately represents th2 geomagnetic noise at 


frequencies of intsrest in MAD operations. 


Be. RECOMMENDATIONS 


A tentative index was developed, tested, and sample data 
obtained. The data analysis for the preliminary svstem was 
accomplished using a mainframe computer. It is reconmended 
“ieectienes: work include improving the index and «he 
setting up of a system to on-line décode the inco 
and utiiizing a desk too computer such as th: 
enable real time determination of the MAD noise i: 

The spatial coherence of MAD geomagnetic noise should 
al 


Gameocila=ion. Meg ittenal ly tae t=esebility of the »priot 


@pbe investigated, with possible application to neoiss 


7) 


estimaticn of geomagnetic noise should be evaluated, 
BCSszbiy by using the prediction aethods ororesed by 
Peaser-Smeth fRef. 42, 43). 
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APPENDIX A 


MAD INDEX DATA ANALYSTS SOFTWARE 
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